Model of molecular bonding based on the Bohr-Sommerfeld picture of atoms 
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We develop a model of molecular binding based on the Bohr-Sommerfeld description of atoms 
together with a constraint taken from conventional quantum mechanics. The model can describe 
the binding energy curves of H2, H3 and other molecules with striking accuracy. Our approach 
treats electrons as point particles with positions determined by extrema of an algebraic energy 
function. Our constrained model provides a physically appealing, accurate description of multi- 
electron chemical bonds. 



<n : 
3: 

Or 

a 

^ ■ 

q ; 
o ■ 

C/3 . 
>v 

43 ■ 



> ■ 
\o ■ 
00 : 

in ■ 
o ■ 
\o : 
o . 

o : 

>>: 
43 : 

Oh- 



X 



Quantum chemistry has achieved excellent agreement 
between theory and experiment by use of computa- 
tional power to provide an adequate description electron- 
electron interactions pj. The conventional treatment 
of molecular structure are based on solving the many- 
particle Schrodinger equation with varying degree of so- 
phistication, ranging from Diffusion Monte Carlo meth- 
ods, coupled cluster expansion, configuration interac- 
tions, to density functional theory All are intensely nu- 
merical, limited to rather small systems and at the ex- 
pense of providing a simple physical picture of the chemi- 
cal bond. Despite the successes of modern computational 
chemistry, there remains a need for understanding elec- 
tron correlations in some relatively simple way so that we 
may describe ground and excited states of large systems 
with reasonable accuracy. 

Our goal here is to advance an intuitively appealing 
model of molecular bonding capable of producing binding 
energy curves at chemical accuracy of a few milli-Hartree. 
Our approach is based on the recently resurrected Bohr's 
1913 model for molecules 2], which is derivable from an 
infinite dimensional reduction of the Schrodinger equa- 
tion [3j. The resulting electron configurations are remi- 
niscent of the Lewis electron-dot structure introduced in 
1916 The surprising feature of our work is that all 
molecular binding energy curves studied below can be ac- 
counted for by mostly electrostatic interaction of "well- 
placed" electrons, as if all the complicated kinetic and 
overlapping integrals have been approximated by their 
mean-values via well-chosen electron positions. Such an 
approach can potentially describe the structural elements 
of large molecules beyond the current capability of ab ini- 
tial methods. 

We will first derive our model for H2, followed by ap- 
plications to more complex molecules. Fig. ^ displays 
various electron distances in H2. Distances and ener- 
gies are in units of the Bohr radius ao (— h 2 /me 2 ) and 
Hartree (= eV^o) respectively. The original molecular 
Bohr model [2j quantize the electron's angular momen- 
tum about the molecular axis resulting in the ground 
state energy function 0, Q , 
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FIG. 1: Electronic distances in H2 molecule. The nuclei A 
and B are fixed a distance R apart. 



where the first term is the Bohr kinetic energy and V is 
the Coulomb potential energy given in terms of electron 
distances defined in Fig. ^ 
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R is the internuclear separation. In our model, electron 
configurations of a physical state correspond to extrema 
of an energy function, such as Eq. 3j,||J. 

In Fig. |21 (curve 1) we plot the ground state E{R) of 
Eq. ifTL together with "exact" quantum mechanical re- 
sults Q (solid circles). The original Bohr model yields 
a fairly accurate H2 ground state energy E(R) at small 
R, but becomes increasely inaccurate at larger internu- 
clear separations. This can be remedied by the following 
observation. At large R each electron in H2 feels only 
the nearest nuclear charge, resulting in two weakly in- 
teracting neutral H atoms. Therefore, at large R each 
electron's angular momentum ought to be quantized rel- 
ative to the nearest nucleus, rather than to the molecular 
axis. This asymptotic consideration yields the following 
H2 energy 
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For R > 2.77 this energy function has a local minimum 
for the top configuration of Fig. |3| We plot the corre- 
sponding E(R) in Fig. 0( curve 2). However, at R < 2.77 
the local minimum of the energy function @ disappears 
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FIG. 2: Potential energy curve of the ground state of the H2 
molecule obtained from the Bohr model with molecular axis 
quantization (curve 1) and quantization relative to the nearest 
nucleus (curve 2). Solid circles are the "exact" energies 0. 



and each electron can collapse onto the other nucleus, 
i.e., ri,i and/or r a2 can vanish. As one can see from Fig. 
13 the energy function @, which is a natural generaliza- 
tion of Bohr's hydrogen atom to the molecular case, is 
in good quantitative agreement with the "exact" energy 
over the range of R > 2.77 where the local minimum ex- 
ists. This encourages us to seek a way of extending the 
applicability of Eq. ® to the entire range of R. 

In the above naive generalization of Bohr's atom to the 
molecular case, each electron can collapse onto the other 
nucleus because there is no corresponding Bohr kinetic 
energy term about that nucleus to prevent the collapse. 
By incorporating further insights from quantum mechan- 
ics, we can remove this instability by a simple constraint. 
Quantum mechanically, the two electrons are described 
by a wave function ^(ri, Y2). Electron 1 is a charge cloud 
with a most probable radius r. Let 
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be the quantum mechanical potential between the elec- 
tron cloud centered at nucleus A and the nuclear charge 
of B, or vice versa. In the Bohr picture we treat the 
electron as a point particle on a sphere of radius r cen- 
tered about nucleus A. A subset of the spherical surface, 
a "circle" of positions r satisfying 



— = 9(r,R) 
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will give the correct quantum mechanical interaction en- 
ergy with nucleus B. Thus if we impose the above as a 
constraint, and choose the electron location only from 



this subset of the positions, then r^i can never be zero, 
because the expectation value in Eq. is finite. 

One can derive the effective potential $(r, R) from any 
simple two-electron wave functions, such as the Heitler- 
London (HL) || or the Hund-Mulliken (HM) ;9f wave 
function 'J'. The HL wave function is 

# = 0(1)6(2) ±6(l)o(2), (6) 

where "+/-" corresponds to singlet/triplet state and 
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for i = l — 2 are variational wave functions with param- 
eter a. If we take a(l) as a variational wave function for 
an isolated hydrogen atom, then the variational energy 
is E = a 2 / 2 — a. This reduces to the Bohr model energy 
function for the hydrogen atom if we identify a = 1/r, 
where r is the radial distance from the nucleus. We will 
also use r to denote the radial distance of an electron 
from its nearest nucleus. 

For the HL wave function Eq. J2J yields the well known 
Coulomb and exchange integral, 
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with overlap S = J a(l)6(l)dri. More explicitly, the sin- 
glet and triplet potentials are respectively 
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with f(r,R) = l/R-exp(-2R/r)(l/r+l/R),g(r,R) = 
exp(— R/r)(l + R/r)/r and 
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The singlet state Hund-Mulliken wave function \& = 
[o(l) + 6(1)] [o(2) + 6(2)] yields the HM singlet effective 
potential: 
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[f(r,R)+g(r,R)}. (11) 



For the triplet state, the HL and HM wave functions are 
the same, yielding identical triplet potentials. Note that 
by introducing singlet and triplet potentials, we have aug- 
mented the original atomic Bohr model with elements of 
Pauli's exclusion principle. The latter is essential for any 
successful description of atomic and molecular systems. 

When we apply the constrained Bohr model to H2, the 
resulting energy function is 
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FIG. 3: Electron configuration for the ground (singlet) and 
triplet states of H2 molecule. 

The energy function has an extremum when r al = — 
r and r a i = r&i. The resulting electron configurations 
corresponding to the ground and triplet states are as 
shown in Fig. [31 where, for singlet ground state 

and for the triplet excited state, 



i?$ 2 (r,i?) R' 

These are just geometric distances between the two elec- 
trons expressed in terms of R, r and r^i . The angle 9 is 
determined by the relation l/$ 2 = R 2 + r 2 — 2rRcos9. 

The binding energy curves E(R) for both singlet and 
triplet states are shown in Fig. 0] There are no fit- 
ting parameters in our calculations. The solid and dot- 
ted lines are results from using the HM and HL poten- 
tial respectively. Solid circles are "exact" results 0. 
The constrained Bohr model gives a surprisingly accurate 
E(R) at all R, yielding a ground state binding energy of 
E B = 4.50 eV for the HL potential and E B = 4.99 eV for 
the HM potential [l(J. The "exact" result is E B = 4.745 
eV pj ■ The Heitler-London- Wang effective charge calcu- 
lation (dashed curves) give s substantially worse accuracy 
with E B = 3.78 eV 0nl y 

more elaborate varia- 
tional calculations with configuration interaction can pro- 
duce energies comparable to our constrained Bohr model 
results. 

Generalizing the constraint Eq. © to a system of 
several hydrogen atoms is straightforward. Let's consider 
electron 1 belonging to its nearest nucleus 1 and denote 
the distances from electron 1 to nuclei i as r% (i = 1, 2, 
. . . ). Then the constraint equation reads 

-E^ = E^i^), as) 

where Ri is the separation between nucleus 1 and nucleus 
i. Mutual spin orientation of electrons 1 and electron i 
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FIG. 4: Potential energy curves of the ground and first 
triplet state 3 E+ of the H2 molecule. Solid lines are obtained 
from the constrained Bohr model with HM effective potential, 
while the small dot line is derived with HL potential. Dashed 
curves are from HL effective charge variational treatment. 

(belonging to nucleus i) determines a singlet or triplet $i 
in Eq. Ijl3|l . In this way, we have incorporated elements 
of Pauli's exclusion principle into our model. 

When applying our model to the triatomic H3 
molecule, we consider linear and triangular configura- 
tions as shown in the insert of Figs. [31 and [S] The spacing 
between the nearest nuclei is assumed to be the same, 
equal to R. Due to symmetry, the central electron in the 
linear H3 molecule must be at equal distances from the 
two neighboring nuclei. For this electron, since its posi- 
tion is fixed, there is no collapse and therefore no need 
for any constraint. We only need to constrain the two 
outermost electrons (see insert of Fig. [SJ via Eq. H13|) in 
the form 

----=$ s (r 1 ,R) + $ t (n,2R), (14) 
ri r 3 

where r-j are defined in the insert of Fig. [5] For the linear 
H3 ground state, adjacent electrons in the molecule have 
opposite spins, requiring the singlet potential <& s (ri, i?). 
(We use the HM singlet potential given by Eq. (fTTJl ). In 
this case, the spins of the two outermost electrons must 
be parallel requiring the triplet potential $ t (ri,2i?.). 
Minimizing the resulting energy function 

E=^ + ^a+V (15) 

yields the solid line potential energy curve of Fig. \5\ 
which essentially goes through the "exact" solid circle 
results. 
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constraint Eq. (|13fl reads 
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FIG. 5: Electron configuration and the ground state E(R) of 
the linear H3 molecule obtained from the constrained Bohr 
model (solid curve) and "exact" numerical solution of the 
Schrodinger equation (solild circles). 
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FIG. 6: Electron configuration and the ground state E(R) of 
the triangular H3 molecule. Solid curve is the result of the 
constrained Bohr model while solid circles are the "exact" 
numerical answer. 



Insert of Fig. shows the electrons' positions for the 
equilateral triangle nuclei configuration. We assume that 
electron 1 has spin opposite to those of electrons 2 and 
3. Symmetry dictates that electron 1 lies above, while 
electrons 2 and 3 lie below the nuclear plane along the 
bisector of the equilateral triangle. For electron 1 the 



while for electrons 2 and 3 we have 
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$ s (r 3 ,i?) + $ t (r 3 ,.R). 



Minimization of the energy function 
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with the constraints ljl6|l and (|17|l results in the solid line 
potential energy curve as shown in Fig. |SJ Again for $ s 
we take the HM effective potential given by Eq. Qllfl. 
As in the case of the linear H3 molecule, the constrained 
Bohr model yielded very accurate E(R) over the entire 
range of internuclear separation. The constrained Bohr 
model also gives good results for other molecules, e.g., 
Be2 as shown in Fig. [7|and H 4 . 

In the vicinity of the energy minimum the constrained 
Bohr model for Be2 provides accuracy of a few milli 
Hartree with no fitting parameters. However, since Be2 
is very weakly bound the binding energy is off by about 
50%. Nevetheless, the bond length remains quite accu- 
rate. 
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FIG. 7: Configuration of outer electrons and the ground state 
E(R) of the Be2 molecule obtained in the constrained Bohr 
model (solid curve) and the "exact" result (solid circles). The 
Bohr model provides accuracy of 0.003 Hartree. 

In summary, we have shown that the atomic Bohr 
model, when supplemented by potential constraints ob- 
tained from quantum mechanics incorporating Pauli's ex- 
clusion principle, gives a remarkably accurate description 
of diatomic and triatomic molecules. This constrained 
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Bohr model provides a physically appealing geometric 
picture of how multi-electron bonds are formed and holds 
promise for future applications to complex many-electron 
systems. Possible application of the model includes the 
simulation of biological molecules, where there is no ab 
initio way of doing the calculations with any other meth- 
ods. The model can also be applied to the calculation of 



potential curves of molecules in a super strong magnetic 
field on surfaces of white dwarf and neutron stars. 
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